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Abstract 

Nb/Cu superconducting superlattices have been fabricated by dc magnetron sputtering. This 
system shows a vortex glass transition with critical exponents similar to high temperatures super- 
conductors exponents. The transition dymensionality is governed by the superconducting coupling 
regime. The vortex glass transition shows a pure two dimensional behavior in decoupled superlat- 
tices and a quasi-two dimensional behavior in the superlattice coupling regime. 
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Since the discovery of High-Tc Superconductors (HTCS) vortex matter Physics has called 
the attention of researchers in many fields, for example the vortex state provides a perfect 
realm to study properties of liquid, crystalline and glassy phases. A plethora of different 
phases is observed in HTCS,^ induced by the interplay of vortex- vortex interaction, thermal 
fluctuations, different kinds of disorder, anisotropy, dimensional effects, etc.^ But the most 
remarkable characteristic of the phase diagram is the existence of two different states: a 
magnetically irreversible zero-resistance state, and a reversible state with dissipative trans- 
port properties. In the absence of disorder, the former corresponds to a vortex-solid phase 
with topological order, and the latter to a vortex-liquid phase, both of them separated by 
a first-order melting transition. In the presence of strong quenched disorder, however, the 
topological order of the vortex-lattice is lost, and the zero-resistance state corresponds to a 
vortex-glass ( VG) , and the transition into the dissipative liquid state becomes a continuous 
second-order phase transition .^'^ 

After the rich phenomenology of vortex- phases came out with HTSC, Low-Tc Super- 
conductors (LTCS) have been seldom revisited to check the possible application of HTCS 
paradigms. The existence of a melting transition in Nb single crystals has been reported, 
but in the case of Nb thin films the existence of a glass transition remains controversial.^'^ 

In this paper, we report on the observation of VG transitions in sputtered Nb/Cu low- 
temperature superconducting super lattices. This system has been chosen mainly because 
the presence of quenched disorder is similar to the one existing in sputtered Nb thin films 
and because the artificially layered structure allows tailoring the anisotropy of the system 
in a controllable fashion. In particular, the coupling of Nb layers through the Cu spacer can 
be easily tuned. This has allowed us studying the dimensionality of the VG transition in 
several different coupling regimes. 

The VG transition has been investigated by measuring the electrical transport properties 
in the mixed state of super lattices, in particular isothermal I-V characteristics in applied 
magnetic fields. These characteristics have been collapsed, according to the scahng rules 
proposed in the VG theory, in terms of critical exponents.^ Besides, the consistency of the 
scaling analysis has been checked with some independent criterion, as recently proposed. 
Moreover, we have found a dimensional crossover for a quasi-2D into a pure 2D VG transi- 
tion, governed by the coupling of Nb layers in the superlattice. 

Nb/Cu superlattices were grown on Si (100) substrates using dc magnetron sputtering at 
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room temperature in Ar atmosphere. Several series of superlattices Cud[Nbd'/Cud]N were 
grown with N being the number of bilayers. Structural characterization was made by X-Ray 
diffraction (XRD) . The structural properties of Nb/Cu superlattices have been investigated 
early by other authors. In our multilayers XRD shows that Cu layers are oriented (111), 
while Nb ones are (110). We have refined the spectra using SUPREX program.^^ From 
refinements the modulation length A, as well as several sources of disorder at the interfaces, 
like roughness or interdiffusion have been obtained. The samples studied here did not present 
interdiffusion, and have moderate roughness at the interfaces, ranging from 0.2 to 0.6 nm, 
being larger when the Nb layers are thicker. The samples were lithographed by wet etching 
into a measuring bridge 1 mm long and 100 /xm wide for magnctotransport experiments with 
standard four-probe configuration. 

Magnctotransport experiments were made in a liquid He cryostat provided with a super- 
conducting solenoid. The superconducting coherence length ^5 as a function of temperature 
was calculated from measured upper critical fields Hc2{T), obtained from magnetoresistance 
measurements at constant fixed temperatures R{H)t- The in-plane coherence length 
(parallel to Nb/Cu interfaces) and the perpendicular one ^^j, have been calculated by using 
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respectively. 
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e^ii (T) = [0o/27r/7e2± {T)r' and ^s± (T) = [(j)oH,2± (T) /27r {H^^h (T)) 
As an example, the observed behavior for sample Nbi3nm/Cu27nm is shown in Fig. 1. The 
perpendicular critical field displays typical linear dependence on temperature Hc2±(T)(x(l- 
T/Tf.) (see Fig. 1 inset). However, the parallel critical field shows up a crossover from linear 
dependence Hc2\\(T)(x(l-T/Tc) at high enough temperatures to square-root Hc2\\(T)(x(l- 
T/TcYf"^ at lower temperatures. A dimensional crossover takes place when the perpendic- 
ular coherence length ^5_l(T) reaches a value of the order of dc^ (thickness of Cu in the 
super lattice). Below the crossover temperature, T2D ~0.5Tc, the Nb layers are decoupled, 
showing up two-dimensional (2D) behavior. ^° In other cases, as for instance Nb3.4nm/Cu2.4nm., 
^si.{T) > dcu at all temperatures, and thus superlattices are always in the coupled regime. 

We have measured I-V characteristics with magnetic field H applied perpendicular to 
Nb/Cu layers. For each fixed value of apphed field H, we measured a set (~20) of isothermal 
1-V curves at decreasing temperatures, as those shown in Fig. 2 and Fig. 3. For all measured 
samples and applied magnetic fields the results are similar. The isotherm at the highest 
temperature displays linear behavior at all current levels, with Ohmic resistance R — V/I^ 
Rn, at this temperature H — i?c2±- For characteristics at slightly lower temperatures, the 
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Ohmic response is observed only up to a threshold current level Ini, above which the curves 
are non- linear. The Ohmic resistance in the low-current limit lim V /I ^ is smaller as the 
temperature decreases, as well as the onset of non-linear response shifts to lower current 
levels as temperature is reduced. Below a given temperature, isotherms become highly non- 
linear within the whole experimental window, and isotherms show up negative curvature in 
the low-current limit yielding zero resistance lim V/I — 0. 

The phenomenology described above suggests the existence of a continuous transition 
from a truly superconducting phase with zero-resistance (VG) to a vortex-liquid dissipa- 
tive phase closer to Hc2- As proposed by Fisher-Fisher- Huse,^ and shown experimentally 
for many HCTS systems, this glass-transition is a second order transition and the physical 
quantities must scale with the VG correlation length ^vG and the characteristic relaxation 
time r. These two magnitudes diverge as temperature approaches the glass transition tem- 
perature Tg, following ^vG oc (T — Tg)^'^ and r oc (T — Tg)"^" , with z and v the dynamic and 
static critical exponents. Scaling laws have been proposed to collapse onto a single master 
curve all I-V (or E — J) isotherms within the critical region, by means of the relation^ 

E^VGT ^ J^^g'C± (JMvG'/kBT) [Eq. 1] 
where D is the dimensionality of the glass transition, (± is a universal scaling function above 
(C+ ) or below (C_ ) Tg, 0o the flux quantum, and ks the Boltzmann constant. We have 
applied this scaling analysis to I-V characteristics measured within the experimental window 
10-^V<V<10-^V and RIO"^ A. The voltage cut-off is 10"^ V, since above this hmit all 
characteristics deviate towards Ohmic behavior. Therefore, the scahng analysis is not longer 
valid. ^-"^'^^'^^ This can be seen in the inset of Fig. 2 (b), where the slopes of isotherms [as 
plotted in Fig 2 (a)], d{logV)/d{logI), are displayed. 

For Nb/Cu superlattices whose Nb layers are coupled, the scaling analysis yields very 
good collapses, as the one shown in Fig. 2 (b). Gritical exponents were around z fa 4 and 
v ~1.8, and the dimensionality was always D—2. These parameters are not magnetic field 
or sample dependent. The values of the critical exponents z and u are within the range 
4-6 predicted by theory.^ Moreover, we want to underline that the obtained values are very 
similar to that observed for quasi-2D VG transitions in HTGS.^^'^'' For instance, the same 
values of the critical exponents were found in YBa2Gu307_5 thin films. This fact stresses 
the universality of the VG transition, earlier suggested. ^^'^^ 

Recently, Strachan et al.^^ have proved that the scaling method used is misleading, even in 
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the case of universal critical exponents with adequate values were found. They have argued 
that experimental limitations related to voltage sensitivity floor might allow achieving good 
collapses with arbitrary values of some of the scaling parameters, as for instance Tg. Thus, 
Strachan et al}^ have proposed that a new criterion to unambiguously determine Tg should 
be met: isotherms above and below Tg, but with equal \{T — Tg) /Tg\, must have opposite 
concavities at the same applied current level. In our Nb/Cu superlattices this criterion was 
checked for every set of I-V curves that had been successfully scaled. In the inset of Fig 2 
(b), the derivatives d {logV)/d (logl) of the I-V isotherms shown in Fig. 2 (b) are plotted. 
From the scahng procedure we obtained Tg=3.994 K for this set of isotherms. As can be seen 
in that Figure, the isotherms above and below this temperature (marked with an arrow in 
Fig. 2(d)) met the criterion. At the lowest current levels, upward and downward isotherms 
are observed at similar distance \{T — Tg)/Tg\ below and above Tg. Furthermore the VG 
transitions theory tells that the isotherm at Tg fulfills the relation E oc J"+^ in the low 
current limit, where a — {z + 2 — D)/{D — 1) } Therefore, using the values obtained in our 
scahng, D—2 and z ^A, the expected slope of the critical isotherm would be o; + 1 5. 
Although the critical isotherm is not in the set of measured characteristics, we can give lower 
and upper limits to its slope from those of the isotherms just above and just below Tg. In 
the low current limit (Fig. 2 (b) inset), the first isotherm below Tg has a slope larger than 
~5, whilst the first one above Tg has a slope smaller than ~4.5 . Thus the critical isotherm 
in between them should have a slope around ~4.5-5, as expected from the independently 
obtained scaling parameters. 

The behavior of superlattices whose Nb layers were decoupled was investigated too. The 
inset of Fig. 3 (a) shows I-V characteristics of sample Nbi3nin/Cu27nm in applied magnetic 
field fioH=0.4 T, such that all isotherms, within the critical region, were at temperatures 
well below T=T2d ~0.5Tc (see Fig. 1). Attempts to use the scahng rule (Eq. 1) did not 
yield collapses as the ones obtained in the coupled regime. In fact, from the derivatives of 
I-V curves [Fig. 3 (b)], one can observe some features that rule out the scaling of these 
isotherms according to a quasi-2D or 3D VG transition. On one hand, it is not possible 
determining a finite Tg using the criterion outlined above, since two isotherms of opposite 
concavities at the same current level can not be found. Besides, as can be seen in Fig. 3 (b), 
a maximun slope a maximum slope of ~9 is displayed by the lowest temperature isotherm 
among those showing up a decreasing slope in the low current limit (marked with an arrow) . 
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The isotherm at Tg should be above this one, and therefore it should have a slope larger 
that ~9. As we said before, in a quasi-2T) or 3D VG transition the slope of the critical 
isotherm is a + 1 {z -\- 2 — D) / [D — 1) + 1, which would imply that the critical exponent 
z>8 for D=2, or z>17 if D=3. These are very high values of the critical exponents, not 
supported by theory.^ As argued earlier, this strongly suggests that a quasi-2D or 3D 
VG transition has to be dismissed. However, a good scaling of the isotherms is achieved by 
assuming a different pure 2D VG transition, as first proposed by Dekker et al}^, and later 
found in HTSC systems. ^^'^^ In a pure 2D VG transition, the glass transition temperature 
Tg—0. Therefore a pure 2D VG phase does not exist at any finite temperature, although 
in-plane correlations (2D) develop, diverging as ^yq oc \jT^' when temperature approaches 
Tg — 0. In this transition, the scaling of the isotherms is achieved by plotting /9ea;p[(To/T)*'] 
vs. J /T^^^' , where p — E/ J is the resistivity, Tq is a characteristic temperature, and p and 
1/^—2 are characteristic exponents of the pure 2D VG transition. The exponent p is related 
to the mechanism of vortex motion: p >1 for thermal activation over the relevant energy 
barriers, whereas p ~0.7 is expected in the case of quantum tunneling across them.^^ As can 
be seen in Fig. 3 (a), a good collapse has been obtained with parameters To=300±20 K, 
p=1.05±0.02, and i/'=2. 

Finally, in our experimental situation, a layered superconductor with magnetic field ap- 
plied perpendicular to layers, one may distinguish between the in-plane VG correlation length 
^VG\\, and the perpendicular one ^vG±{^^ong the vortex line).^ The quasi-2D character of the 
glass transition was explained in high Tc superconductors assuming that anisotropy induces 
limited vortex length, that precludes ^vg± to diverge .^^'^^ Thus, when approaching Tg, only 
^VG\\ diverges up to the macroscopic size of the sample, whereas ^vg± would remain finite 
with nearly a constant value. This applies to Nb/Cu superlattices in the coupled regime. 
Following Yamasaki et al}^ and Zefrioui et al}"^ we can estimate an upper limit of ^vGi. from 
I-V characteristics. Isotherms above Tg show up Ohmic behavior at low current level, but 
they becomes non-linear above Ini- At this current level, the work done by the Lorentz force 
to create vortex excitations equals the thermal energy ,^ Jnifj^oivGwivGL = ksT. We used 
the isotherm at the highest temperature within the critical region, in particular the one at 
T=4.038 K shown in Fig. 3 (a), for which !^125 Acm~^. The in-plane correlation length 
^VG\\ should be larger than the mean inter-vortex distance Oq = (0o/A*o-f^)^^^ ~ 130 nm for 
/iqH — 0.11 T. In particular we assumed ^vg\\>'^o,q, and thus we obtained ^yG±<90 nm. 
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That is, the correlation length along vortex line is always shorter than sample thickness. 
However, it may be longer than other relevant characteristic lengths, as the superlattice 
modulation length A==40 nm and the superconducting coherence length Cs±(T)^ 35 nm. 
We have estimated the vortex length / in the regime where this superlattice is decoupled, 
in which we observed a pure 2D VG transition, with the work done by Lorentz force 
Jni(l>o^vG\\l = ksT}^ Taking the isotherm at r=2.150 K [Fig. 3 (a)], J„/ ^80 Acm'^ , and 
^yG||>2ao =150 nm for iiqH — 0.4 T, we get I <30 nm. Therefore, the vortex length I is 
shorter than superlattice modulation length A=40 nm, and cannot be much longer than the 
coherence length at this T, ~ 20 nm. A picture emerges from those estimations, in 
which coherence length ,^5^ and sample thickness d are the relevant length scales to which 
correlation length along vortex line ^VG± has to be compared. The quasi-2D character of the 
glass transition in the coupled regime develops since CvG±, longer that $,s±, is shorter than 
sample thickness d. At lower temperatures, Nb layers are decoupled by Cu ones. Therefore 
the vortex length I is hmited below modulation length A. Because of this, the vortex length 
/ and coherence length ^5_L.are similar This yields a pure 2D VG transition. This is similar 
to what is observed in highly anisotropic HTCS, for which pure 2D VG transitions have 
been observed when vortex length is limited to superconducting coherence length. 

In summary, we have shown strong evidence of the VG transition in the mixed state of 
low-temperature superconducting Nb/Cu super lattices. The HTCS analysis apphes directly 
to this system, and the same universality has been observed, in spite of the very small thermal 
fluctuations of LTCS in comparison with HTCS. Besides, a dimensional crossover from a 
quasi-2D into a pure 2D VG transition has been observed, which is governed by the ratio of 
the VG correlation length ^vg to the superconducting coherence length ^g. 
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Figure 1 



FIG. 1: Superconducting coherence lengths of sample Cu27nm[Nbi3nm/Cu27nm]ioas a function on 
temperature, both parallel ^^[[(T) and perpendicular '^5_l(T) to Nb/Cu layers. Inset: Measured 
parallel (black circles) and perpendicular (open circles) critical fields Hc2 • Solid lines are linear fits 
Hc2±(T)(x(l-T/Tc), while the dashed one is the best fit to Hc2\\(T)(x(l-T/Tcy/^. 
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Figure 2 
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FIG. 2: (a) I-V isotherms for sample Cu27nm[Nbi3nm/Cu27nm]io in applied field jLtoH=0.11 T at 
temperatures (from left to right) 4.060 K>T>3.865 K, separated '^5-15 mK. The dashed line 
separates isotherms above and below Tg=3.994 K. Vertical and horizontal dotted lines delimit 
the experimental window used for scaling, (b) Scaling of the above isotherms as explained in the 
text . Inset: Derivatives of the log(I)-log(V) isotherms at temperatures (from bottom to top) 
4.060 K>T>3.908 K. Black dots are within the experimental window used for scaling. The arrow 
separates derivatives of isotherms just below and above Tg. 
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Figure 3 
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FIG. 3: (a) Scaling of the isotherms shown in the inset with a pure 2D VG model. Inset: I-V 
isotherms for sample Cu27nm[Nbi3nm/Cu27nm]io in applied field /ioH=0.4 T at temperatures (from 
left to right) 2.160 K>T>1.780 K, separated ~5-15 mK.. Horizontal dotted line delimit the ex- 
perimental window used for scaling, (b) Derivatives of the log(I)-log(V) isotherms at temperatures 
(from bottom to top) 2.160 K>T>1.856 K. Black dots are within the experimental window used for 
scaling. The arrow marks derivatives of the isotherm at the lowest temperature showing decreasing 
slope in the low current limit 
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